Nontrapping arrest of Langmuir wave damping near the threshold amplitude 
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Evolution of a Langmuir wave is studied numerically for finite amplitudes slightly above the 
threshold which separates damping from nondamping cases. Arrest of linear damping is found to 
be a second-order effect due to ballistic evolution of perturbations, resonant power transfer between 
field and particles, and organization of phase space into a positive slope for the average distribution 
function f av around the resonant wave phase speed Vj,. Near the threshold trapping in the wave 
potential does not arrest damping or saturate the subsequent growth phase. 
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Plasma theory has usually been pursued independently 
of the theory of critical phenomena. Recently, how- 
ever, it has been revealed that evolution of a monochro- 
matic electrostatic Langmuir wave of finite amplitude in 
a Maxwellian plasma is a threshold phenomenon. Specif- 
ically, after a short initial period of approximately linear 
damping according to Landau's classic theory, a wave 
with initial amplitude Aq greater than a threshold 
stops decreasing and starts to grow approximately expo- 
nentiallybefore undergoing irregular oscillations in am- 
plitude 0,0] (Fig.Q. Both the amplitudes and times at 
which the wave first ceases to damp and grow (labelled 
"arrest" and "saturation" ) are power-law functions of the 
difference (A - A*) 0, thus casting the process into a 
new universality class of dynamic critical phenomena. 

For a collisionless plasma the distribution function 
(DF) is usually not Gaussian, and because of the long- 
range character of the Coulomb force these systems are 
outside the domain of equilibrium thermodynamics. Un- 
like the theory of critical phenomena in thermodynamics, 
where only spatial correlations are considered through 
the order parameter 0(x) and the partition function 
Z = J T>(j>{x)e~ H ^ 4], threshold physics in collisionless 
systems involves correlations in velocity space || . There- 
fore these correlations, due to resonant energy exchange 
between particles and waves, are a new paradigm for crit- 
ical phenomena potentially applicable in a vast class of 
systems, e.g. coupled phase oscillators which show Lan- 
dau damping or equivalent Josephson-j unction arrays Q. 

Crawford's pioneering analysis reveals the striking 
difference between thermodynamic and plasma situations 
due to this physics: the resonance between particles and 
waves at the phase velocity vj, = w pe / k turns the thermo- 
dynamic exponent /3 = 1/2 |4| into the "trapping scaling" 
exponent (3 = 2, which describes saturation of the weak 
bump-on-tail || and gravitational instabilities]^ . In the 
frameworks of linear and quasilinear theory [9j arrest of 
the linear damping of plasma waves (as well as saturation 
of the growth) might be explained in terms of flattening 
of the DF at v<j,, thus bringing the damping (growth) rate 
7l ~ (df /dv) v=V(j> of a kinetic instability to zero. 

Trapping of electrons in a monochromatic wave's elec- 



tric potential is often suggested as a nonlinear mecha- 
nism to stop the initial exponential damping phase and 
to saturate the wave's growth 0,0]. Trapping and its 
associated Bernstein- Greene-Kruskal (BGK) modes 0] 
also imply a certain shape of the DF plus trapped and 
untrapped orbits in velocity phase space. However, it is 
controversial whether trapping is relevant to the damp- 
ing threshold. For instance, one analysis 0] assumes 
ergodicity of trapped particles in a single-wave potential 
and predicts the threshold initial electric field amplitude 
Eq through the critical ratio q c — |7L|/^tr ~ 0.06 of 
the absolute Landau damping rate |tl| to the trapping 
frequency mtr = {kEoe/nie) 1 / 2 . In contrast full Vlasov- 
Poisson (V-P) simulations for a Maxwellian plasma yield 
q c w 0.85 from the asymptotic evolution and q c ~ 1.0 
from the initial evolution 0, with constants of propor- 
tionality slightly different from unity for other thermal 
plasmas 0- 

Other conflicting evidence exists on the role of trap- 
ping. Consider the critical exponents T mm , /3 m i n , r sa t and 
/?sat for the power-law functions of (A — Aq) obeyed by, 
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FIG. f: (a) Two-component electron (in blue) and ion (red 
fine) m = 1 field evolution for m p /m e — 1836. (b) Electric 
field envelope for mode m = 1 for immobile ions (black solid 
line) and mobile ions (blue dashed line), and the ratio of the 
simulated frequency to its analytic prediction, ui* /lu (right 
axis, black dash-dotted line). Diamonds mark the "arrest" 
time r, mm and "saturation" time t sa ^. 
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respectively, the time t m - m and amplitude A m i n at which 
the initial damping phase finishes, as well as the time 
i sa t and amplitude A S3i ^ at which the first exponential 
growth phase saturates y|: e.g., £ m j n oc (A — j4Q)~ r min 
and A m i n cx (A — A^ram. First, the temporal expo- 
nents r min = 0.901 ± 0.008 and r sat = 1.039 ± 0.011 are 
measurably different from each other and the value 0.5 
expected from the definition of u>t T . Second, the field ex- 
ponents /3 sat = 1.88 ±0.07 and f3 min = 2.72 ±0.09 are re- 
markably different from each other and the value /?tr = 1 
expected for trapping [llj. These points argue against 
trapping causing either of the arrest and saturation phe- 
nomena. Third, calculations with Ao ^> Aq lead to r and 
f3 exponents closer to 0.6 and 1.3, respectively, and the 
oscillation spectrum has clear peaks near wtr, suggesting 
that trapping plays a role well above threshold |3j. 

In this Letter we first simulate one-dimensional (1- 
D) V-P two-component plasma with initially Maxwellian 
distributions for electrons and ions and demonstrate that 
ion mobility does not affect the threshold phenomenon 
for Langmuir wave damping seen in V-P simulations 
without ions. Then, using one-component electron V- 
P simulations, we demonstrate that the DF phase por- 
trait when the wave first ceases to damp is much simpler 
than a BGK equilibrium [l2T | and shows no evidence for 
trapping. Instead, we demonstrate that the initial DF 
resonantly evolves a positive slope in velocity space that 
stops the initial Landau damping and supports the sub- 
sequent exponential growth. We also demonstrate that 
the DFs are different at the arrest and saturation times 
and are not consistent with trapping. 

To clarify the importance of ion mobility we employ 
first the two-component 1-D V-P model, normalizing to 
electron quantities: 

df a /dt + vdf a /dx-ii a Edf a /dv = , (1) 

/+oo 
(f p -f e )dv. (2) 
-oc 

Here a = e,p, m a is the particle mass, f a the compo- 
nent's DF, [i e = I, [ip = —m e /m p , and E{x,t) is the 
electric field. The boundary conditions are assumed to 
be periodic. The initial electron distribution is 

f e (x,v,0) = l/v / 27rw the exp(-w 2 /2^ he )[l+A cos(fc m a;)] , 

where v t ^ e is the Maxwellian thermal speed for electrons, 
Ao the initial electric amplitude, k rn — 2-KrajL is the 
wave number of the mode m, and L is the length of the 
system. The ions are initially uniform and Maxwellian- 
distributed in velocity space with T p = T e . 

The simulations use m = 1, «the — 0.4, Debye length 
A De ~ 0.31, and L = 2tt ps 20.18A De - They have 
N x = 256 cells in the x direction both for electrons 
and ions, and N ve = 20000 and N vi = 2000 cells in 
speed for electrons and ions, respectively, within the 
domains [— 10i>th a i lO^tha]- The Cheng-Knorr method 



[ijj was used to solve Eqs Q and J5J) with double 
precision. System invariants I^ a = J fadxdv are con- 
served better than A/ 3e //3 e (0) < 10~ 6 for electrons, and 
AI 3p /I 3p (0) < 10" 9 for ions. 

Fig- Q] shows the evolution of the mode m = 1 for ini- 
tial amplitude A Q = 0.012, A*, = (8.51 ± 0.06) x 10" 3 , 
and m p /m e = 1836. This type of evolution is observed 
experimentally |l6j |. The existence of significant ion mo- 
tion in Fig. Ufa) seems, at first glance, to suggest that 
the evolution is seriously affected by ion mobility. How- 
ever, the envelope field amplitude of the electron oscil- 
lations in Fig. ^a) is almost identical to that for immo- 
bile ions [Fig.^b)]. Quantitatively, the initial damping 
phase in Fig. ^ stops at time t m - m « 4410^5 and ampli- 
tude A m i n ~ 1.64 x 10 -5 , and is then followed by almost 
exponential growth which saturates at t sa ± ~ 1365 lo~^ 
and A ss x w 2.42 x 10~ 4 . These quantities are identical to 
those calculated in the electron V-P simulations of Ref. 
[3j, where m = 4 was assumed for the perturbation and 
u th = 0.1 for the electron thermal speed. This is ex- 
pected because fcADe, the wave frequency u>, and 7l are 
the same for the two simulations. 

Analytic theory predicts that to ~ 1.2851 a-> po , but 
the simulated value uj* — 1.2705 ± 9 x 10~ 4 is slightly 
shifted from w due to the large value of Ao and varies 
slightly with time [Fig. Gib)]. Linear damping rate is 
7l i=a — 0.0661 uj pe . For smaller Aq — 10 -5 both uj* and 
7l match the standard Landau theory very well (not 
shown), with {\(u*- u)/u\, |( 7 £ - 7l)/7l|} < 2 x lO" 4 . 

These two-component V-P results demonstrate that 
the threshold phenomenon for Langmuir wave damping 
is robust against ion effects. Accordingly one-component 
simulations, with ions acting as a neutralizing back- 
ground, are used below. 

The DF near the phase velocity = uj* jk\ « 1.271 
at these moments is shown in Fig. |21 and reveals dras- 
tic discrepancies between the evolution which ends with 
arrest of damping at t = £ m i n , and the subsequent evo- 
lution until the growth saturates at t — i sa t. At the 
moment t = i m j n the phase space portrait reveals no 
signs of particle trapping - only filamentation due to 
phase mixing (Fig. the upper view). Moreover, in- 
stead of a stationary state this distribution supports ap- 
proximately linear (meaning exponential) growth on the 
interval i m ; n < t < i sa t, as Fig. ^b) shows. Crucially, 
the DF at £ sa t does not consist of the closed orbits (or 
whorls in velocity-position space) expected for trapping. 
Instead, the orbits are still open, although they clearly 
indicate progress towards trapping. Trapping is there- 
fore responsible for neither the arrest of damping nor the 
saturation of the growth phase. 

In the linear theory developed by Landau 1] growth is 
due to a positive slope in the DF at the phase velocity 
of the wave, \v\ — v^. Fig. shows the DF averaged 
on x coordinate, fo(v,t) = (l/£) J Q f{x,v,t) dx at t = 
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FIG. 2: DF when (upper panel) t = t m ; n and damping stops 
and (bottom panel) t = t 8a ^ and growth saturates. 



i m ; n . Instead of the flattening of fo near the resonant 
velocities v — ztv^ predicted by quasilinear theory 0, 
/o(w,£ m in) acquires a positive slope in a small vicinity 
of 1)0, and therefore can support (approximately) linear 
growth after the moment t = t m - m as Fig. ^b) shows. 

Contrary to the situation near t — £ m i n when damping 
ceases and the physics looks quite smooth and regular, 
fo becomes quite irregular near the time t = t sa t when 
growth saturates (Figs.[5]and0J). In particular, the lower 
panel of Fig. is strongly reminiscent of trapping, al- 
though strictly closed trajectories do not appear for this 
Aq. Also, while on average the slope of /o(u, f sa t) at 




FIG. 3: fo(v, t m i n ) (solid line) and fo{v,0) (dashed line) 
on semilogarithmic (upper panel) and linear (bottom panel) 
scales for two velocity intervals: (a) —2.0 < v < 2.0 and (b) 
1.24 < v < 1.32. 




FIG. 4: fo(v,t sa t) (solid line) and fo{v, 0) (dashed line) 
on semilogarithmic (upper panel) and linear (bottom panel) 
scales for two velocity intervals: (a) —2.0 < v < 2.0 and (b) 
1.24 < v < 1.32. 



v = dzVtj, seems to have decreased compared with time 
t = ^min [Fig- EI a )]? if varies irregularly in the neigh- 
borhood of ±i>0 and therefore may support excitation of 
oscillations with a wide range of phase speeds. 

Fig. shows the evolution of, and power transfers 
between, the average DF fo(v,t) and the DF compo- 
nents f\(v,t) and /2(u,i) at k\ and k2, respectively, with 
\f m (v,t)\ = {Re 2 [f m {v,t)\ +lm 2 [f m {v,t)}} 1 / 2 . It shows 
that the dynamical picture can be divided into regions 




FIG. 5: Evolution of DF components fo (upper panel), /i 
(middle panel) and fa (bottom panel). Red lines are DF com- 
ponents at the moments t m ^ n and t sa i, e.g. \fi(v, i sa t)|. 
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FIG. 6: Power transfer P(v) for a damping wave with 7 = 0.03 
in Eq. (15) of Ref. 0). 



with distinct characteristics that identify the processes 
causing the evolution. Fig. shows that the turbulent 
processes responsible for the (relative) flattening of /o in 
the resonant area near v<f> start only after t SSu t , when spa- 
tial Fourier components E m other than m = 1 become 
comparable to E\ (not shown here). 

The ripples of /o, |/i|, and I/2I in time and velocity 
appear to be "fingerprints" of ballistic change of ini- 
tial perturbation and power transfer between the field 
and particles. The latter claim is justified by Fig. El 
which illustrates the power transfer rate for a wave 
growing/damping linearly by resonant wave-particle in- 
teractions given by Eq. (15) in Ref. [17|. Together 
with Fig. 13 which shows the evolution of Sfo(v,t) — 
[fo(v,t) - fo(v,0)}/f (v,0) on the interval < t < t min , 
Figs El and Fig. clearly demonstrate that the physical 
process responsible for arrest of linear damping is the res- 
onant power transfer between the wave and the m = 
and higher order components of the DF. 

An insight into the striking difference between the crit- 
ical exponents (3 m i n and /3 sa t comes from critical phenom- 
ena theory: critical exponents depend on the properties 
of correlations for a specific system (e.g., on its dimen- 
sionality) and/or a universality class (e.g., Ising, perco- 
lation, surface growth etc.) [4j. The DFs in full phase 
space (position and velocity) are different at times i m i n 
and i sa t (see Figs 2-5 and Fig. EJ), so the critical ex- 
ponents might be different. This difference is contrary 




to the idea that trapping explains both the arrest and 
saturation phases, which should result in the same ex- 
ponents. Some plausibility for velocity-space structures 
having this effect follows from 1-D V-P self-gravitating 
calculations: varying the resolution in v seriously affected 
estimates of the "trapping scaling" exponent (3 — 2 5] . 

In summary, we studied the V-P model for initial 
Langmuir wave amplitudes slightly above the thresh- 
old that separates damping and non-damping evolution. 
Electron-ion simulations show that ion mobility does not 
modify the threshold found for Langmuir damping in 
electron-only simulations. Phase space diagnostics show 
no signs of trapping or the DF flattening near t — ^min 
instead the combined effects of ballistic evolution of per- 
turbations and resonant power transfer at |i>| ~ are re- 
sponsible for arrest of the linear (Landau) damping then. 
Since the spatially-averaged DF is not flat at t m [ n but in- 
stead has a positive slope near the resonant velocity v$, 
this state is not stationary but instead leads to (linear) 
growth which is saturated at t = t 8a t. The saturation 
time i S at marks the boundary between the regular and 
stochastic evolution of the wave electric field, again with 
no evidence for trapping saturating the growth phase. 
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FIG. 7: Evolution of 5fo{v,t) till the moment t = t T 
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